Particle rearrangements during transitions between local minima of the potential 
energy landscape of a supercooled Lennard-Jones liquid 
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The potential energy landscape (PEL) of supercooled binary Lennard-Jones (BLJ) mixtures ex- 
hibits local minima, or inherent structures (IS), which are organized into meta-basins (MB). We 
study the particle rearrangements related to transitions between both successive IS and successive 
MB for a small 80:20 BLJ system near the mode-coupling temperature Tmct- The analysis includes 
the displacements of individual particles, the localization of the rearrangements and the relevance 
of string-like motion. We find that the particle rearrangements during IS and MB transitions do not 
change significantly at Tmct- In particular, an onset of single particle hopping on the length scale 
of the inter-particle distance is not observed. Further, it is demonstrated that IS and MB dynamics 
are spatially heterogeneous and facilitated by string-like motion. To investigate the mechanism of 
string-like motion, we follow the particle rearrangements during suitable sequences of IS transitions. 
We find that most strings observed after a series of transitions do not move coherently during a 
single transition, but subunits of different sizes are active at different times. Several findings suggest 
that, though string-like motion is of comparable relevance when the system explores a MB and when 
it moves from one MB to another, the occurrence of a successful string enables the system to exit a 
MB. Moreover, we show that the particle rearrangements during two consecutive MB transitions are 
basically uncorrelated. Specifically, different groups of particles are highly mobile during subsequent 
MB transitions. We further find the positions of strings during successive MB transitions weakly 
but positively correlated supporting the idea of dynamic facilitation. Finally, the relation between 
the features of the PEL and the relaxation processes in supercooled liquids is discussed. 

PACS numbers: 64.70 



I. INTRODUCTION 



For supercooled liquids approaching the glass transi- 
tion temperature Tg, the viscosity increases continuously 
by more than 13 orders of magnitude. Despite valuable 
progress in recent years, a complete theory that ratio- 
nalizes the corresponding slowing down in molecular dy- 
namics is still lacking. On the one hand, mode coupling 
theory (MCT) predicts a power-law divergence of relax- 
ation times at a critical temperature Tmct^- Though 
such a temperature dependence provides a satisfactory 
description of experimental data above Tmct'^j a sin- 
gularity at the critical temperature is not observed and 
Tmct is typically 1.1 — 1.5 Tg. This deviation is usually 
attributed to hopping processes that restore ergodicity 
below Tmct, but are not included in the ideal MCT. 



On the other hand, potential energy landscape (PEL) 
approaches, based on the pioneering work of Goldstein^, 
have proven useful in the field of the glass transition phe- 
nomenon^*^. Here, the high-dimensional vector of all 
particle coordinates is considered as a point moving on 
the potential energy surface. A PEL description is par- 
ticularly useful at sufficiently low temperatures where, as 
proposed by Goldstein^, the evolution of the system can 
be decomposed into vibrations about local minima of the 
PEL on short timescales and transitions among distinct 
minima on long timescales. Then, knowing the proper- 
ties of the local minima of the PEL is often sufficient to 
calculate observables. 

In computational work on supercooled liquids, a de- 
tailed picture of the PEL is accessible. Following Gold- 
stein's idea, Stillinger and Weber^ introduced the concept 
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of basins where a basin in configuration space is defined 
as the set of points that via steepest descent path along 
the PEL maps onto the same local minimum, or inherent 
structure (IS). The resulting tiling of the configuration 
space into non-overlapping basins of attraction simpli- 
fies the thermodynamic description of the system. For 
example, it is possible to write the free energy approxi- 
mately as a function of the energies and the vibrational 
frequencies of the IS^ so that the statistical properties of 
these minima determine all thermodynamic quantities. A 
link between thermodynamics and kinetics is suggested 
by the Adam-Gibbs relation which connects the relax- 
ation times with the configurational entropjjiS. Such a 
link was qualitatively established by Sastry et alpii who 
observed that the onset of non-exponential relaxation is 
accompanied by the sampling of IS with progressively 
lower potential energy upon cooling. Later, the Adam- 
Gibbs equation was shown to hold quantitatively in sim- 
ulations of watev^, a binary Lennard-Jones (BLJ) liq- 
uid^^. and silicai^. Further, a dependence of the fragility 
on the shape of the PEL was reportedi^ii^. Such a rela- 
tion was long proposed by Stillinger who suggested that 
strong liquids have a uniformly rough PEL, while the IS 
of fragile liquids are organized into meta-basins (MB)M, 
or, equivalently, funnelsi^*^. 

Several workers studied the relation between the prop- 
erties of the PEL and the trajectories of the individ- 
ual particles in more detail. Schr0der et alm^ confirmed 
Goldstein's picture^ showing that the dynamics of a BLJ 
liquid below T « Tmct can be separated into local vi- 
brations and transitions between IS. Further, string-like 
displacements of groups of particles were found during IS 
transitions similar to the motion observed in the equilib- 
rium liquidi^. Finally, single-particle hopping indicated 
by a secondary peak in the van Hove correlation function 
at T < Tmct was reported not to result from the crossing 
of single energy barriers, but from sequences of IS tran- 
sitions''-^. Instantaneous normal mode analysis suggests 
that a transition from non-activated to activated dynam- 
ics takes place at Tmct^^'^^'^^'^^- At T > Tmct, 
the system is always close to a multi-dimensional ridge 
between two basins, and the slowing down of dynam- 
ics arises from the reduction of directions along which 
free exploration of the configuration space is possible. 
In contrast, at T < Tmct, the system samples regions 
of the configuration space that have no free directions 
so that activated processes are required for relaxation. 
However, the method to determine a correct distribution 
of saddle points is still a matter of debate^^^. In dif- 
ferent approaches, no indication for a change of the dy- 
namical behavior at Tmct was observed, but activated 
barrier crossing was found to be relevant already above 
Tmct^^'^^^- Doliwa and Heue n^^i^^ calculated the dif- 
fusion constant D{T) from the static properties of the 
PEL by taking into account the organization of the IS 
into MB. In this way, the effect of back-and-forth jumps, 
which frequently occur between the IS within the same 
MB, but do not contribute to long-range diffusion, can be 



eliminated, and only hopping between MB is considered, 
which should resemble a random walk on the PEL. It 
was further shown that a sequence of jumps is required 
to escape from low-energy MB. According to Denny et 
al.^'' "transitions between IS within a MB involve very 
small flexing of a cage, while transitions between MB in- 
volve collective rearrangements". Wales and coworkers 
distinguished between "non-diffusive" and "diffusive" re- 
arrangements depending on whether nearest neighbors 
are changed and found that the energy barriers involved 
in both types of rearrangements differ more in strong 
than in fragile liquids^^. Further, they observed that 
upon cooling the distance atoms move during an IS tran- 
sition decreases significantly at T kTmct^- 

Despite this progress, the relation between the dynam- 
ics of single particles in real space and the motion of 
the system in configuration space is still elusive. In par- 
ticular, the gap separating MCT and PEL approaches 
has not yet been bridged^. In this article, we present 
a detailed study of the particle rearrangements result- 
ing from transitions between local minima of the PEL 
of a supercooled BLJ liquid. Due to the organization of 
IS into MB, we consider transitions between consecutive 
IS and consecutive MB where the latter result from a 
series of the former. For simplicity, we refer to the corre- 
sponding displacements as IS dynamics and MB dynam- 
ics, respectively. Stillinger^ related MB dynamics to the 
a-process and IS dynamics to the (Johari-Goldstein) (3- 
process'^^ and, hence, our approach may yield valuable 
insights into the nature of these relaxation processes in 
supercooled liquids. The present analysis of IS/MB dy- 
namics includes measures for the displacements of the in- 
dividual particles, the localization of the rearrangements 
and the relevance of string-like motion, which is believed 
to facilitate the structural relaxation in supercooled liq- 
uidsiSiSiiSSi. Wc focus on two temperatures Th > Tmct 
and Ti <Tmct- In this way, we can investigate whether 
a change of the transport mechanism observed for the 
equilibrium liquid at T « Tmci "i^'^i^^i^^^ cf. section |lll 
manifests itself in distinct IS/MB dynamics. Further, by 
following the particle displacements during suitably de- 
fined sequences of IS transitions, we study the mechanism 
of string-like motion in detail. To analyze IS /MB dynam- 
ics beyond the concept of single transitions, we investi- 
gate the correlation of the particle displacements during 
consecutive IS/MB transitions. In view of our results and 
recent experimental findings, we suggest a modification 
of Stillinger's picture of the relation between the features 
of the PEL and the relaxation processes in supercooled 
liquids. 

II. MODEL AND SIMULATION 

The simulations used to generate the data analyzed 
here are described in ReffS^. In summary, we investigate 
a BLJ liquid characterized by the interaction potential 

Vo,p{r) = Asa^p [{cT^ti/ry^ - [ao^p/rf] (1) 



with the parameters N = Na + Nb = 65, Na = 52, aAB = 
0.8 a AA, cfbb =0-S>ScrAA, £ab = ^-^ £aa, £bb = 0-5 Eaa 
and = 1.8. Linear functions are added to the potentials 
to ensure continuous forces and energies at the cutoff re- 
These modifications of the original potential of Kob and 
Andersen25i^ are necessary for the simulation of small 
system32fliS£. The data were generated using Langevin 
molecular dynamics simulations (MD) with fixed step 
size A = 0.015 = (2A:_BrAi/mc)^/^, equal particle masses 
TO, friction constant ^ = 1 and periodic boundary condi- 
tions. All results are reported in units of aAA, "m-, £aa 
and nv^X^ /2eAA for length, mass, energy and time, re- 
spectively. 

At regularly spaced times t;, we quench the MD tra- 
jectory x{t) to the bottom of the visited basins;^-, yield- 
ing the discontinuous trajectory of IS, In the IS 
trajectory, the entries at subsequent ti are identical un- 
til a transition to a new basin of attraction takes place. 
These multiple entries are eliminated so that the resulting 
trajectory only contains the IS separated by transitions. 
Further, we apply an interval bisection method described 
in Refj2^ to ensure that all relevant transitions between 
the regular quenches are resolved. The final set of trajec- 
tories contains all successively visited IS and con- 
sists of more than 10^ IS for the studied temperatures 
T/ = 0.435 and T,, = 0.500. From ^{tj), we extract the N 
trajectories of the individual particles, f^(tj), which form 
the basis of the present analysis. In addition, a straight- 
forward algorithm is applied to construct the MB^SiSi. 
First, all time intervals [ijjij] are searched where is 
the time of the first and i^- ^ tj the time of the last oc- 
currence of the IS S.{tj). Next, any two overlapping time 
intervals [tj,tj] U are combined. A similar proce- 
dure was used in Ref>2^ to calculate the diffusion constant 
D{T) from the properties of the constructed MB. Strictly 
speaking, IS that occur exactly once in the trajectory 
(,{tj) {tj =t'j) form trivial MB. Here, we exclusively use 
the term MB for non-trivial MB. Thus, the MB tran- 
sitions studied take place between non-trivial MB and 
consist of several IS transitions. 

For PEL approaches, it is essential to use small sys- 
tems, because otherwise interesting effects are averaged 
Q^^j37j38j39^ On the other hand, the system should not be 
too small so as to avoid significant finite size effects. As 
was discussed previouslyS*^, = 65 turns out to be a 
good compromise. Specifically, comparing the dynamics 
of systems with A^ = 65, A^ = 135 and N = 1000 it was 
observed that the respective diffusion constants D{T) are 
identical within 20% for all temperatures T > Tmct^- 
Further, a critical temperature rMCT = 0.45±0.01 is ob- 
tained from a power-law fit to D{T) for A^ = 65 as com- 
pared to a value Ta/ct = 0.435 estabhshed for larger BLJ 
system3Si2&. Finally, the behavior of the A^= 130 system 
resembles that of two independent copies of the A^ = 65 
system^. Therefore, significant finite size effects are ab- 
sent for T « TmcTi and - similar to what is known for 
larger BLJ systems^^*^ - the A^ = 65 BLJ system can be 




0.2 0.4 0.6 0.8 1 

FIG. 1; Probability distributions p (dris) and p {drMB) char- 
acterizing the displacements of the A particles during transi- 
tions between consecutive inherent structures and consecutive 
meta-basins, respectively. Data for 0.435 and = 0.500 
are compared. Lines: Exponential decays obtained by inter- 
polating the data at drjs/MB^^-'^- 



regarded as a model of a typical supercooled liquid. 

In previous work on the studied system^Si^, Doliwa 
and Heuer calculated the mean waiting time in the MB, 
< tmb [T) > ■ Its relation to the time constant of the a- 
process, Tq., depends on the system size. For A^ = 65, a 
ratio Ta/ <tm_b>~ 30 is found at T; and T^. Further, 
the mean waiting times in the IS, <t/s(T)>, are ap- 
proximately a factor of six shorter than the correspond- 
ing < tmb{T) >. For the diffusion constants, a ratio 
D{Th)/D(Ti)^ 10 was reported ^5.26^ However, for vari- 
ous models of equilibrium liquids, the transport mecha- 
nism changes significantly in a comparable temperature 
range. Specifically, a secondary peak develops in the van 
Hove correlation function upon cooling T — > Tmct , which 
is usually interpreted as an onset of single-particle hop- 
pjjjgiiji7j33j34^ Moreover, Th lies well in the tempera- 
ture regime where the idealized version of the MCT is 
valid^Si^, whereas a strong contribution of activated pro- 
cesses proposed by the extended version of this theory 
can be expected at T/ < Tmct- Therefore, despite a 
moderate variation of the diffusion constant, the nature 
of the structural relaxation may change significantly in 
the studied temperature range. 



III. RESULTS 
A. Transitions between inherent structures 

First, we study the particle displacements drjs result- 
ing from IS transitions, i.e., the displacement of parti- 
cle i during the transition ^{tj) — > ^(ij+i) is given by 
dris = \'''^{tj+i) — r"^{tj)\- In Fig. 1, the probability dis- 
tribution p (dris) for the A particles of the BLJ liquid is 




FIG. 2: Probability distributions p{dris;dEis) describing 
the displacements of the A particles during transitions be- 
tween consecutive inherent structures with energy difference 
dEis (r = 0.435). We distinguish between energy differences 
dEis < 0.5 (dEio), 0.5 < dEis < 1.5 (d^me) and 1.5 < dEis 
(dEhi). Lines: Exponential decays Ai exp{—dris/dro) where 
the factors Ai were chosen to match the respective curves 
p {drjs; dEjs) and drg was obtained from an interpolation of 
p{dris), cf. Fig. 1. 



displayed. For both temperatures, a rapid decay domi- 
nates the distributions. Specifically, for sufficiently large 
displacements, p{dris) decays exponentially. Thus, the 
distributions for T > Tmct and T < Tmct do not ex- 
hibit a striking difference, but the temperature depen- 
dence manifests itself in variations of the mean particle 
displacement (<(ir/s(ri)> = 0.081, <drisiTh)> = 0.093). 
These findings are consistent with those of Schr0der et 
alii on a 50:50 BL J mixture {N = 800) where p {drjs) for 

— 5/2 

T « Tmct was described by a power-law drjg at small 
dris and by an exponential decay at large dris- There, 
it was concluded that the latter functional form results 
from particles taking part in a local event, whereas the 
former is caused by particles adjusting to this event. Con- 
sistent with this interpretation, the deviations from an 
exponential behavior at small drjs are less pronounced 
in our case of a smaller system. In contrast to Schr0der 
et al^, we - except for an analysis of string-like motion 
- do not include the B particles, because the dynamics 
of the minority component of the 80:20 mixture, though 
qualitatively similar, is somewhat fasteriSS*^ so that a 
superposition of the respective results may disturb the 
underlying functional forms of the distribution functions. 

The IS transitions are further characterized by the dis- 
placement of the whole configuration, drfg = Y!i^=i dials' 
and by the absolute value of the energy difference of 
the involved IS, dEjs- For T; and T^, we find that 
the probability distribution p {dEis) is well described by 
an exponential decay at sufficiently large dEis- More- 
over, the mean energy difference decreases upon cooling 
{<dEisiTi)> =0.78, <dE I siTh)> =0.98). The probabil- 



FIG. 3; Probability distributions p {dris', Eis) characterizing 
the displacement of the A particles during transitions start- 
ing from an inherent structure of energy Eis at T = 0.435. 
We differentiate between energies Eis < —297 (Eio), — 297< 
Eis < -295 (E^i), -295 < Eis < -293 (Emh) and -293 < Eis 
{Eh,). 



ity distributions p {drfg) for the studied temperatures are 
broad and asymmetric. They exhibit a peak at drfg «4 
and extend to drfg « 17, indicating that the IS dynamics 
shows a large diversity. Thus, an analysis in terms of 
statistical quantities is necessary. 

Next, the relation between the IS dynamics and the 
energy difference dEis is analyzed by means of the prob- 
ability distributions p {dris', dEis). In Fig- 2, we see that 
the mean displacement of the A particles increases with 
dEis. For Ti, < dris> = 0.06 and < dris> = 0.13 are 
obtained for energy differences dEis < 0.5 {dEio) and 
dE> 1.5 {dEhi), respectively. Thus, a large variation of 
the potential energy is accompanied by large particle re- 
arrangements. To study IS dynamics in different regions 
of the PEL separately, we consider the probability dis- 
tributions p {dris; Eis) where Eis is the energy of the 
initial IS of the transition. The results for Ti in Fig. 
3 show that the energy Eis weakly affects the particle 
displacements. Hence, IS dynamics in different regions 
of the PEL are comparable. We note that the results 
for Tfi are qualitatively similar to those shown in Figs. 2 
and 3. Further, it was checked that the basic features of 
our findings do not depend on the specific choice of the 
energy limits. 

It is known from both experiments^is^Sii^ and simu- 
lations^i*22ii4i^ that molecular dynamics in supercooled 
liquids are spatially heterogeneous. We analyze whether 
particles experiencing similar displacements due to IS 
transitions are spatially correlated, too, as suggested by 
Schr0der et ali^. Since a cluster analysis is limited by the 
system size N = 65, we focus on the A particles with the 
largest and the smallest dris during the respective tran- 
sition and analyze the mobility of their neighbors. Specif- 
ically, we first rank all particles according to their dris 




FIG. 4: Probability distributions pirnis) and p{mMB) for 
the mobilities of the particles in the next neighbor shells 
of the A particles with the largest (mo) and the smallest 
(im) displacement during the respective inherent structure 
and meta-basin transition. See text for details. Dashed line: 
p(Tn/s/Mfl) = l/64. 




FIG. 5: Typical example of string-like motion during a se- 
quence of transitions between inherent structures. The gray 
and the black spheres visualize the particles before and af- 
ter the sequence, respectively. The different shades of gray 
denote particles replacing their neighbors at different times. 
It is seen that the "macro-string" after the series contains 
two "micro-strings" due to individual transitions. The micro- 
string on the left hand side is active first, the one on the right 
hand side 5 inherent structure transitions later. 



and assign mobilities mis so that m/5 = 65 is attributed 
to the most mobile and mis = 1 to the most immobile 
particle. Then, we pick the most mobile A particle, re- 
arrange where required the mobihties of the remaining 
particles (to/s = 1,...,64) and calculate the probabil- 
ity p {mis) that a particle in the first neighbor shell of 
the selected particle has the mobility mis- Finally, this 
procedure is repeated for the most immobile A particle. 
Hence, if mobile and immobile particles were randomly 
distributed throughout the system a probability distri- 
bution p{mis) — 1/64 would result for both the most 
mobile and the most immobile A particle. Fig. 4 shows 
that a random distribution does not apply. Instead, the 
particles with the largest displacements are mostly sur- 
rounded by other highly mobile particles. Thus, spatially 
heterogeneous dynamics is also observed when analyzing 
IS dynamics of a 65-particle BLJ liquid close to Tmct- 

For supercooled liquids close to Tmct, it has been 
demonstrated that the relaxation of highly mobile par- 
ticles is facilitated by "string-like motion" which means 
that groups of particles follow each other in quasi-one- 
dimcnsional pathsiSiSiSiiSiil, see Fig. 5. In addition, 
such strings were reported for the IS dynamics of a 
50:50 BLJ liquidii. Here, we quantify the contribution 
of this dynamical pattern to the IS dynamics in more 
detail. Similar to Donati et sXwiS,^ we construct strings 
by connecting any two particles i and k if min[ — 
r''{tj+i)\,\r'-{tj+i) -r'^{tj)\] < 5 = 0.6, where we now 
consider both the A and the B particles of the mixture. 
Since 5 is smaller than the hard-core radii of the A and 
the B particles, this condition implies that one particle 
has moved and another particle has occupied its position. 
With this definition, string-like motion is observed dur- 
ing 29% of the transitions at T/. On the other hand, for 



86% of those transitions where string-like motion occurs, 
the most mobile particle is involved in a string. This 
is consistent with a spatial correlation of highly mobile 
particles, cf. Fig. 4. 

The strings can be further characterized by their length 
Z, i.e., the number of participating particles. In Fig. 6, we 
display the mean number of strings of length I involved 
in one transition, < nis{l)>. This quantity is related 
with the average number of particles moving in strings 
during one IS transition by < TV/s > = J2i "-/slO >• 
For Ti and Th, the decrease of < nis{l)> is consistent 
with an exponential decay. However, we find that the 
functional form of the decay depends somewhat on the 
choice of S. Moreover, finite size effects can be expected 
for I > 5 due to the small system size. In any case, ex- 
ponential decays were also observed for the probability 
distribution of the string length when analyzing the dy- 
namics of equilibrium liquidsiSi^i*^Siii. We emphasize 
that the observed temperature dependence of < TV/s > 
does not imply that string-like motion is more important 
at higher temperatures. Instead, one must take into ac- 
count that only particles that show a certain minimum 
dris can fulfill the above specified replacement criterion 
and, thus, a larger mean displacement - as was observed 
at the higher temperature, cf. Fig. 1 - tends to lead to 
a larger number <A/7s>. Indeed, it was found for equi- 
librium liquids that the fraction of particles moving in 
strings increases upon coolingiS*^. 

Finally, we study the localization of particle rearrange- 
ments resulting from IS transitions. For this purpose, we 
measure the number of involved particles by the quanti- 
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FIG. 6: Mean number of strings of length I observed during 
transitions between successive inherent structures (< n/s(Z)>) 
and meta-basins {<nMB{l)>), respectively. Data for temper- 
atures Ti =0.435 and 0.500 are shown. Lines: Interpola- 
tions with an exponential decay. 



ties zijs and Z2./5. The former is defined as 
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where dRjs is the maximum displacement of an A par- 
ticle during the respective transition, and the latter is 
calculated according to2S4& 
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In the case that n particles move the same distance and 
the remainder is immobile, zijs and 2:2,7s equal n. In 
Fig. 7, we show the probability distribution p (zijs) for 
the studied temperatures. For T; and Th, the distribu- 
tions are nearly symmetric and peak at zijs ~ 16. The 
shape of all curves can be described by a Gaussian with 
a width parameter a « 4.5, which again reflects the di- 
versity of IS dynamics. The distributions p (22, /s) (not 
shown) are close to a Gaussian centered at Z2js ~ 17 and 
characterized by a~7.5. A closer inspection reveals that 
the mean values <zijs> and <Z2js> decrease by about 
0.5 when decreasing the temperature from Th to Tj, i.e., 
IS dynamics becomes slightly more local upon cooling. 
We emphasize that the value of about 16 for the number 
of particles involved in IS dynamics should not be taken 
too literally since measures of the localization different 
from zijs and Z2js can yield values that are up to a 
factor of 3 smaller. 
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FIG. 7: Probability distributions p(zi,/s) and p(zi,mb) 
characterizing the number of particles participating in the 
rearrangements due to transitions between successive inher- 
ent structures and meta-basins, respectively, cf. Eq. |5| The 
results for T; = 0.435 and Th = 0.500 are compared. Dashed 
lines: Interpolations of the data for T with a Gaussian. 



B. Transitions between meta-basins 

Recalling that Stillingeri^ related the a-process to MB 
transitions, we now turn to the particle rearrangements 
during these transitions. For this purpose, we search the 
IS with the lowest energy in each MB, ^(tfc), where k is 
the index of the MB, and define the particle displacement 
during a MB transition as the one obtained by a compar- 
ison of ^{tk) and ^(ifc+i). Since k and A; -1-1 are not the 
indices of consecutive IS, but of consecutive MB, these 
displacements are caused by a sequence of IS transitions. 
In what follows, we describe the particle displacement 
during a MB transition by the vector drMB{k) and de- 
note its absolute value as dvMB- 

The probability distributions p{drMB) calculated for 
Ti and T^ are included in Fig. 1. Since MB dynamics re- 
sults from several IS transitions, the mean displacement 
<drMB{T)> is larger than <drjs{T)>. For the tempera- 
ture dependence of MB dynamics, we find different mean 
values < drMB{Ti)> = 0.18 and < drMB{Th)> = 0.20, 
but the shape of the distributions p (drMs) at T; and Th 
is comparable. Similar to p{dris), the distributions are 
well described by an exponential decay at large drujB- 
Interestingly, for neither type of transitions, there is en- 
hanced probability at drjg^f^jg « 1.0. Hence, single- 
particle hopping on the length scale of the inter-particle 
distance observed for various larger models of equilibrium 
liquids at T « Tmct^^^^^'^^^ does not manifest itself in 
IS and MB dynamics of our small system. The mean en- 
ergy difference of the IS involved in the MB transition, 
<dEMB>, increases from 1.42 at Tj to 1.56 at Th- 

To analyze whether particles showing similar mobil- 
ities during MB transitions are spatially correlated we 



calculate the probability distributions p^mMs) which, 
analogous to p{mjs), characterize the mobility of the 
neighbors of the most mobile and the most immobile A 
particle during a MB transition. It is evident from Fig. 4 
that the most mobile particle during a MB transition is 
prevailingly surrounded by other mobile particles. More- 
over, a comparison of p{mis) and p{mMB) shows that 
the spatially heterogeneous nature of IS and MB dynam- 
ics is very similar at T « Tmct- To gain further valu- 
able insights, we study the localization of the particle 
rearrangements due to MB transitions using the quan- 
tity zi^MB defined in analogy to Eq. 12 Inspecting the 
probability distributions p(zi,mb) for T; and Th in Fig. 
7, we see that a comparable number of particles par- 
ticipates in IS and MB dynamics, suggesting that basi- 
cally the same group of particles is "active" during all IS 
transitions involved in the respective MB transition. Fur- 
ther, similar to IS dynamics, MB dynamics is slightly less 
local at the higher temperature {< zi^mb{Ti)> = 17.7, 
<zi^MB{Th)>= 18.3). Hence, we find no indication 
that the length scales attributed to the spatially het- 
erogeneous nature of IS and MB dynamics, respectively, 
increase upon cooling. 

We further explore the spatial heterogeneities associ- 
ated with MB dynamics by investigating the relevance 
of string-like motion. With the same definition as in the 
previous section, we find that string-like motion occurs 
during 87% of all MB transitions at Ti where the aver- 
age number of particles that move in strings during a 
MB transition amounts to <Nmb> =7.2. These num- 
bers indicate that string-like motion yields an important 
contribution to the particle rearrangements during MB 
transitions. In other words, the group of particles that 
takes part in MB dynamics usually achieves the large 
displacements by means of string-like motion. The rel- 
evance of this type of motion for MB dynamics is also 
obvious from Fig. 6 where we display the distributions 
<«a/_b(0> characterizing the mean number of strings of 
length I during one MB transition. 



C. Sequences of inherent structure transitions 

Due to the organization of the IS into MB, two differ- 
ent processes are relevant when the system samples the 
PEL. First, the exploration of a MB, i.e., the back-and- 
forth jumps between the IS in the same MB and, second, 
the transition between distinct MB. In this section, we 
study the particle rearrangements involved in these pro- 
cesses by following the particle rearrangements during 
suitable sequences of IS transitions. Since the diffusion 
constant is basically determined by the trapping of the 
system in long-lived ME^^-^^, we now focus on MB within 
which at least six IS transitions occur, i.e., we consider 
49% of all MB. To investigate the exploration process 
and the transition process separately we construct E- and 
T-sequences using the following criteria: An E-sequence 
combines all IS transitions connecting the two IS within 
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FIG. 8: Probability distributions p^{drsQ) and p'^(drsQ) 
describing the displacements of the A particles during se- 
quences of inherent structure transitions that occur when a 
meta-basins is explored (E) and when a transition between 
different meta-basins takes place (T), respectively (r = 0.435). 
See text for details. For comparison, we included the cor- 
responding results for the involved individual transitions, 
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the same MB that show the largest distance dr'fg. This 
means that, on average, nine IS transitions form an E- 
sequence for T/ . To obtain the T-sequences we unite four 
IS transitions that occur when the system moves to an- 
other MB. More precisely, if two successive MB are sep- 
arated by more than four IS transitions, we choose four 
jumps in the middle of this series. On the other hand, if 
the MB are separated by less than four IS transitions, we 
add the adjoining transitions within the involved MB to 
the sequence. We do not regard the latter procedure as a 
serious problem because, in any case, long-lived MB are 
exited by a series of events^^ so that it is not completely 
clear which individual transitions are a part of the escape 
process. The chosen numbers are also motivated by the 
overall displacements during the E- and T-sequences, see 
below. However, we ensured that the basic features of 
our results are unchanged when the numbers are varied 
in a meaningful range. The IS transitions involved in the 
E- and T-sequences, respectively, are denoted as E- and 
T-transitions in what follows. 

Fig. 8 shows the probability distributions p^'^ (drsg) 
and p^''^ (drjs) for the displacements of the A parti- 
cles during the E/T-sequences and the E/T-transitions 
at T;, respectively. While different mean displacements 
< drjs >^> < drjs >^ are found for the single tran- 
sitions, the rearrangements during the E- and the T- 
sequences are described by nearly identical distribution 
functions. Thus, the smaller mean displacements during 
the E-transitions and the back-and-forth jumps between 
IS within the same MB lead to the effect that four T- 
transitions result in the same overall displacement as, on 
average, nine E-transitions. In other words, the particle 
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FIG. 9; Mean number of strings of length I resulting from: (i) 
single transitions between successive inherent structures vis- 
ited during the exploration (<n/s(Z)>^) and the transition 
process (<n7s(0>^)5 (^0 sequences of transitions of inherent 
structures visited during the exploration (<nsQ(i)>^) and 
the transition process {<nsQ{l)>'^) (r = 0.435). See text for 
details. Lines: Interpolations with an exponential decay. 



rearrangements during the exploration and the transition 
process are very similar at the level of sequences of tran- 
sitions. Despite the quantitative differences, the distri- 
butions for the single transitions, p^{dris) and p"^(c?r/5) 
are still comparable. In particular, due to the correlation 
between the particle displacement and the energy differ- 
ence dEjs, cf. Fig. 2, some of the deviations result from 
the fact that larger mean energy differences are found for 
the T-transitions. Strictly speaking, when calculating the 
distributions p ^ (dr/s ; dEjs) a,nd p'^ {dr i s ; dEjs) we find 
that both distributions are very similar for IS transitions 
characterized by high energy differences, but somewhat 
different in the case of small dEjs- 

The finding {drsQ)^p^ (drsq) enables us to study 
the relevance of string- like motion for the exploration and 
the transition process without any effect of the respective 
mean particle displacement, see the discussion of Fig. 6. 
This was one reason for the specific choice of the param- 
eters when constructing the sequences. The distributions 
<nsQ{l)> ^''^ for the mean number of strings of length I 
during one E- and one T-sequence at Ti , respectively, are 
displayed in Fig. 9. The good agreement indicates that 
string-like motion is of similar importance for the parti- 
cle rearrangements during the exploration and the tran- 
sition process, respectively. This finding excludes that 
it is simply the occurrence of string-like motion that al- 
lows the system to escape from a MB. This point will be 
further discussed later in this section. For comparison, 
we include the distributions < njs{l)> ^''^ characteriz- 
ing string-like motion during the individual E- and T- 
transitions, in Fig. 9. Both distributions can be interpo- 
lated by an exponential decay, suggesting that the basic 



features of string-like motion are comparable for the dif- 
ferent types of transitions. As was discussed above, the 
quantitative differences result at least in part from the 
distinct mean particle displacements due to the E- and 
the T-transitions, cf. Fig. 8. Taking also into account 
that a comparable number of particle participates in E- 
and T-transitions (<2i,/s>^= 16.3, <zijs>^ = 16.3 
for Ti) we conclude that the different types of transitions 
show at most some quantitative differences in the vicinity 
of T]\icT, supporting our prior finding that IS dynamics 
in different regions of the PEL are comparable, see Fig. 
3. 

Next, we study the mechanism of string-like motion. 
When comparing <nsQ{l)>^''^ with their counterparts 
for the single transitions <njs{l)> , cf. Fig. 9, it be- 
comes obvious that the distributions have different slopes 
in a semi-logarithmic representation. Hence, one may 
speculate that a string observed during a sequence does 
not result from a single transition, but from the interplay 
of particle displacements due to subsequent jumps. For 
example, one can imagine that a long "macro-string" , 
defined as a string during a sequence, is composed of 
several "micro-strings" resulting from the displacements 
during the individual transitions of the series. This will 
now be studied in more detail. For this purpose, we con- 
sider all macro-strings of length / > 3 and analyze how 
many of them result from a single transition. We find 
that only 15.4% and 13.7% of the macro-strings during 
the E- and the T-sequences, respectively, are the conse- 
quence of a single event. On the other hand, 48.4% (E) 
and 57.5% (T) of the macro-strings are composed of at 
least one micro-string. The remainder, namely, 36.2% 
(E) and 28.8% (T), can be traced back to concerted 
single-particle type displacements in successive jumps. 
Likewise, for 47.8% (E) and 53.8% (T) of all particles 
involved in macro-strings of length I > 3, the replacement 
of the neighboring particle takes place during a single 
transition, while, for the remainder of the particles, sev- 
eral displacements must add up for the replacement cri- 
terion to be fulfilled. Of course, these numbers depend 
on the definition of the strings. Nevertheless, they show 
that most macro-strings, especially, the long ones, do not 
result from a coherent motion of all particles, but from 
subsequent motions of single particles or small groups of 
particles. 

A typical example of the interplay of micro-strings and 
single-particle motions in forming a large macro-string 
is shown in Fig. 5, where the different shades of gray 
denote particles replacing each other at different times. 
It is also seen that the replacements do not start and end 
at the "head" (right hand site) and the "tail" (left hand 
side) of the string, respectively, but occur in a random 
order. Though ordered sequential replacements along the 
string are observed in many instances, we still regard the 
scenario in Fig. 5 as a typical example. 

Especially, the macro-strings of the E-sequences may 
be a subtle result of individual motions. Since the IS 
within a MB are visited several times, the micro-strings 
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of the E-sequences may show back-and-forth jumps. To 
quantify this effect we calculated the probability that a 
micro-string observed during an E-transition jumps back 
during a later transition in the same MB. We find that 
for 44.5% of the micro-strings all particles jump back 
to their initial positions. Hence, back-and-forth jumps of 
micro-strings are a frequent phenomenon when observing 
IS dynamics within a MB. In comparison, the back-jump 
probability for the micro-strings during the T-sequences 
amounts to only 2%. Considering also the results in Fig. 
9, one may speculate that, though string- like motion is of 
similar relevance for the particle rearrangements during 
the exploration and the transition process, respectively, 
it is the occurrence of successful strings that enables the 
system to escape from a MB. 



D. Correlation of successive particle displacements 

As aforementioned, results of Doliwa and IIeueri2SiS& 
suggest that jumps between MB resemble a random- walk 
on the PEL. To check this conlusion, we study the corre- 
lation of particle displacements resulting from MB tran- 
sitions at two different times, i.e., four-time correlation 
functions are considered. To measure for how many MB 
transitions an A particle remembers the direction of an 
initial motion we define the correlation function 



D{Ak) 



drMB (fc) drMB{k + Ak) 
drMsik) drMB{k + Ak) 



(4) 



Here, the brackets < . . . > denote the average over all MB 
transitions and all A particles. Due to the properties of 
the scalar product, D{Ak) will be positive (negative) if, 
on average, the displacements of a particle during the MB 
transitions £,{tk) and ^{tk+Ak) ^{tk+Ak+i), 

respectively, have the same (opposite) direction. Another 
property of the particle displacement during a MB tran- 
sition is the mobility. Different from the definition used 
so far, we now characterize the relative mobility of a par- 
ticle by 



Kk) = 



drMBjk)- <drMB{k)> 
<drMB{k)> 



(5) 



where < drMsik) > is the mean displacement of an A 
particle during the MB transition ^{tk)^^{tk+i)- Hence, 
the correlation function 



M{Ak) = {fi{k)fi{k+Ak)) 



(6) 



relates the relative mobilities during different MB tran- 
sitions. To study the correlation of the particle displace- 
ments resulting from different IS transitions, we define 
the correlation functions D{Aj) and M{Aj) in analogy 
to their counterparts for the MB transitions. 

Fig. 10 (a) shows ^(Afc) and M(Afc) for the MB dy- 
namics at Ti . Inspecting the latter correlation function a 
rapid decay is evident. For the former, the negative sign 



of D{Ak = — 1) implies that subsequent motions are back- 
ward correlated, i.e., the particles tend to move towards 
their old positions during the next jump. However, the 
values |£>(Afc = l)|<l and Af(Afc = l)<l indicate that 
the particle rearrangements during consecutive MB tran- 
sitions show only a small correlation. Hence, a random 
walk on the PEL is indeed a good approximation for the 
jumps between MB. In addition, the findings imply that, 
during different MB transitions, different groups of par- 
ticles are highly mobile and form strings, corroborating 
our speculation that the occurrence of a successful string 
results in the escape from a MB. 

For comparison, D{Aj) and A/(Aj) characterizing IS 
dynamics at T/ are displayed in Fig. 10(b). Due to the 
back-and-forth jumps between the IS within a MB the 
correlations persist for several IS transitions. In partic- 
ular, sequences of back-and-forth jumps between two IS 
result in a oscillatory behavior of D{Aj). M{Aj) does 
not completely decay within the first 15 IS transitions. 
Since the same transition may occur several times when 
the system explores a MB, the information about the ini- 
tial mobility is not completely lost until the MB is exited. 

Very recently, Garrahan and Chandler^-'-'^ proposed 
that bulk dynamics in supercooled liquids can be un- 
derstood based on two simple ingredients, namely, the 
existence of spatially heterogeneous dynamics and the 
facilitation of dynamics in the vicinity of mobile regions. 
To check this assumption, we study whether strings dur- 
ing consecutive MB transitions are spatially correlated. 
For this purpose, we first characterize the positions of 
all strings by their center of mass. Then, for every string 
during the MB transition k ~> k+1, the minimum distance 
between its position and the positions of any string dur- 
ing the transition A: — 1 — > /c is determined. Finally, we 
repeat this calculation for randomly chosen positions of 
the strings during the MB transition k^k + 1. For the 
actual strings at Ti , we obtain a mean minimum distance 
<d > — 1.3 which is slightly smaller than <d > = 1.5 
found for a random distribution, suggesting that dynam- 
ics in the vicinity of mobile regions may be somewhat fa- 
cilitated. However, the observed effect is weak and needs 
to be validated for larger systems. 



IV. DISCUSSION AND SUMMARY 

We studied the particle rearrangements during transi- 
tions between consecutive inherent structures (IS) and 
consecutive meta-basins (MB) of a supercooled 80:20 
BL J liquid. Specifically, the displacements of the individ- 
ual particles, the localization of the rearrangements and 
the string-like motion were characterized for two temper- 
atures, Th > Tmct and T; < Tmct- In addition, the IS 
dynamics were analyzed in dependence of the energies 
of the involved IS and their energy difference. Consid- 
ering that the motion of the system on the PEL can be 
decomposed into the exploration of MB and the transi- 
tion between distinct MB, we compared the particle rear- 
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FIG. 10: (a): Correlation functions that relate the displacements of single A particles during different transitions between 
meta-basins at T = 0.435. D{Ak) measures the correlation of the directions of subsequent motions; M(Afc) correlates the 
respective relative mobilities, cf. Eqs. 2] and |S| (b): Analogously defined quantities D{Aj) and A/(Aj) characterizing the 
correlation of the particle displacements during different transitions between inherent structures. 



rangements resulting from both processes. This analysis 
was done both at the level of single transitions (E- and 
T-transitions) and at the level of suitable sequences of 
transitions (E- and T-sequences). Finally, time corre- 
lation functions of the displacements during transitions 
between successive MB were analyzed. Since the time 
constants for the dynamics of the A and the B parti- 
cles of the BLJ mixture are somewhat differenfe^S*^, we 
mostly focused on the motion of the former. However, 
with respect to all studied quantities, only minor differ- 
ences are observed for the B particles. 

Comparing the results for >Tmct and T/ <Tmct, 
no change of the basic features except only gradual vari- 
ations were found. For IS and MB dynamics, the mean 
energy differences, <dEjs/MB>, the mean displacement 
of the A particles, < dris/MB>i and the mean number 
of A particles taking part in a transition, < Zijs/mb >i 
decrease upon cooling, but the shape of the correspond- 
ing distributions is basically unchanged. For T/ and T^, 
the distributions p [zijs/mb] are close to Gaussian and, 
over a wide range, the curves p (drjs/MB) decay exponen- 
tially. The decrease of < dEj g /jv/^ > is consistent with 
the decline of fc^T. The variation of < dris/MB> is, at 
least in part, a consequence of the changing mean energy 
difference, because it was found that large changes of the 
energy are accompanied by large particle displacements. 
Hence, when observing IS and MB dynamics in a 65- 
particle BLJ liquid, there is no evidence that the mecha- 
nism for the particle motion is discontinuously altered at 
Tmct- In particular, we do not find single-particle hop- 
ping on the length scale of the inter-particle distance. 

These results are in agreement with findings by 
Schr0der et alii^, who concluded that single-particle hop- 
ping observed for the equilibrium liquid at T ?» Tmct 
does not result from transitions over single barriers, but 
"the jump occurs via a number of intermediate IS" . On 
the other hand, at variance with the outcome of the 



present work, Hernandez-Rojas and Walea^ reported 
a "distinct change in behavior" at T « Tmct when 
investigating the particle rearrangements in a kinetic 
Monte-Carlo approach. Since T/ — 0.435 is still close 
to Tmct — 0.45 ±0.01 it may be useful to check our 
results over a broader temperature range. With this 
reservation, the absence of a significant change of the 
dynamical behavior at Tmct appears to support that 
the sampling of the PEL changes gradually'^'"'-^^ instead 
of discontinuously as concluded in normal-mode analysis 

In addition, we analyzed the dependence of IS dynam- 
ics on the energy of the initial IS, Ejg. Weak variations 
of < drjs{Ejs)> indicate that the energy has a minor 
influence. These results suggest that IS dynamics in dif- 
ferent regions of the PEL are comparable. Such behavior 
is in accordance with the current picture of the PEL of 
fragile glass formers where a comparable ruggedness is 
assumed throughout the landscape^. Moreover, a weak 
dependence of IS dynamics on Eis is consistent with our 
findings for the particle rearrangements during E- and 
T-transitions, i.e., for IS dynamics occurring when a MB 
is explored and when a transition between different MB 
takes place, respectively. We observed that the mean dis- 
placement < drjs> , is somewhat larger than < drjs> ^, 
but this effect results in part from larger mean energy dif- 
ferences during the T-transitions, because generally the 
particle rearrangements are larger for higher energy dif- 
ferences. In addition, the number of particles involved in 
IS dynamics and in string-like motion are similar during 
the exploration and the transition process. Hence, with 
respect to the accompanying particle rearrangements, no 
principal differences exist between IS transitions near the 
bottom of a MB and near the ridges between distinct MB. 

We further demonstrated that the particles partici- 
pating in IS and MB dynamics, respectively, i.e., par- 
ticles showing high displacements are not randomly dis- 
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tributed, but rather reside prevailingly in the next neigh- 
bor shell of each other. Thus, consistent with previous 
computational work on supercooled liquid o^^i^^'^^i'^^i^^ , 
IS and MB dynamics at T w Tmct are spatially hetero- 
geneous. In experiments^i*^, a spatially heterogeneous 
nature was observed for dynamics near Tg. Though the 
dynamical heterogeneities at T ^Tmct and at T>Tg, re- 
spectively, show similar features, e.g., a comparable rate 
memory^^, their exact relation is still elusive. 

The clustering of mobile particles found for IS and MB 
dynamics is in accordance with the occurrence of string- 
like motion. A closer analysis revealed that this dynam- 
ical pattern is observed during 29% of the IS transitions 
and during 87% of the MB transitions at Tj. The latter 
value together with <7Vmb> =7.2 obtained for the mean 
number of particles forming strings per MB transition in- 
dicate that string-like motion is very important for MB 
dynamics. On the other hand, this type of motion ap- 
pears to be less relevant on the shorter time scale of IS 
dynamics. However, one has to consider that during an IS 
transition, the displacement of most particles is too small 
to satisfy the criterion used to define strings. Hence, 
a limited relevance is a natural consequence. To gain 
further insights, we compared string- like motion during 
sequences of IS transitions that take place when the sys- 
tem explores a MB and when it moves from one MB to 
another, respectively. Constructing E- and T-sequences 
such that comparable particle displacements result, we 
observed <nsQ(/)>^« <nsQ{l)>'^ for the mean num- 
ber of strings of length / during the respective sequences. 
Thus, string-like motion during the exploration and the 
transition process is not different, suggesting that the 
probability to find a string depends merely on the dis- 
placements of the particles. 

Concerning the mechanism of string- like motion, it was 
demonstrated that the formation of strings results from 
the concerted interplay of the particle rearrangements re- 
lated to subsequent IS transitions. In other words, most 
macro-strings defined as strings that result from a se- 
ries of transitions do not arise from a coherent motion of 
all particles at the same time, but from coordinated dis- 
placements of subunits at different times. These subunits 
include both single particles and micro-strings, i.e., small 
groups of particles replacing each other within one tran- 
sition. Further, the subsequent particle replacements in 
a macro-string do not necessarily start at the "head" of 
the string and end at the "tail", but successive replace- 
ments often take place at random positions in the string, 
cf. Fig. 5. Moreover, for IS dynamics inside the MB, 
the formation of a macro-string frequently includes back- 
and-forth jumps of the involved micro-strings. Hence, 
all these findings elucidate that string-like motion is a 
complex dynamical process in which several rearrange- 
ments within a group of particles occur. In view of this 
result, the observation <zijs> ~ <zi,mb> for the num- 
ber of particles taking part in IS and MB dynamics, re- 
spectively, becomes plausible. Since particles that are 
"active" during a MB transition mostly participate in 



strings, the multi-step nature of this type of motion has 
the consequence that basically the same group of par- 
ticles is active during all IS transitions involved in the 
respective MB transition, i.e., for several IS transitions, 
high particle mobility persists in a certain region of the 
system. Based on the multi-step nature of string-like 
motion, Gebremichael et &\M^ speculate that the occur- 
rence of strings is related to local structural properties 
of supercooled liquids. For example, they suggest that 
temporary fissures open a quasi-lD channel in which the 
particles can move in strings. Further investigation on 
this point is underway^. 

Finally, we studied the correlation of particle displace- 
ments during different transitions between MB, i.e., time 
correlation functions of displacements were investigated. 
Analyzing both the direction of the motion and the rel- 
ative mobility we find that the displacements of any one 
particular particle during consecutive MB transitions are 
basically uncorrelated. In other words, different groups of 
particles are mobile during successive MB transitions. In 
contrast, when relating the particle rearrangements due 
to different IS transitions, the back-and-forth jumps be- 
tween the IS within a MB and the process of string forma- 
tion result in a correlation of the particle displacements 
that persists for several IS transitions. We also inves- 
tigated whether the strings during two consecutive MB 
transitions are spatially correlated. For this purpose, we 
calculated the minimum distance between a new string 
and any of the old strings. We found a mean minimum 
distance that is 15% smaller than in the case of a ran- 
dom distribution of the new strings, suggesting that new 
strings tend to be formed near regions where string-like 
motion has taken place. This appears to support Garra- 
han and Chandleri^SiS who stress the importance of dy- 
namic facilitation, i.e., they assume that regions showing 
high mobility assist neighboring regions to become mo- 
bile. However, the effects observed here are weak and 
need to be validated for larger systems, as well as in the 
equilibrium liquid. 

As aforementioned, Stillingeri^ assigned the a- and 
the (Johari-Goldstein) /3-process to MB and IS dynam- 
ics, respectively. From experimental work, it is known 
that, different from the a-process, the /3-process exhibits 
an Arrhenius-like temperature dependence^ so that both 
processes often merge into a single high temperature re- 
laxation near Tmct-^- Multi-dimensional NMR experi- 
ments have shown that the molecular reorientations as- 
sociated with the a 54i55,56 ^^^^ /J-processSiSiS close 
to Tg are complex multi-step processes. For the lat- 
ter relaxation, the multi-step reorientation is restricted 
to a small section of the unit sphere^^iSiSiffi where 
the accessible solid angle strongly grows with increas- 
ing temperature^fliSiS. In addition, computational stud- 
ies have revealed that the length scale associated with 
the spatially heterogeneous nature of the structural re- 
laxation in supercooled liquids increases upon cooling 

We now discuss Stillinger's view in the context of the 
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present and the summarized previous findings for super- 
cooled liquids. Several results are inconsistent with the 
assumption that the a-process is identical with MB dy- 
namics. First of all, we found that MB dynamics results 
in a mean particle displacement that is much smaller than 
the inter-particle distance, cf. Fig. 1, indicating that, at 
the studied temperatures, the structural relaxation is not 
complete after a single MB transition. Likewise, Tq is ap- 
proximately a factor of 30 longer than the mean waiting 
time in the MB <tmb>^^^- Finally, the temperature 
dependence of the mean number of particles involved in 
a MB transition, < zi^mb >, is incompatible with an 
increasing length scale attributed to the spatial hetero- 
geneities of MB dynamics. Concerning the /3-process, 
some of the experimental results appear to be at variance 
with the assumption that this relaxation is related to sin- 
gle transitions between IS. Specifically, the experimental 
finding that the /9-process results from a multi-step pro- 
cess is inconsistent with a single jump in the PEL. More- 
over, the strong growth of the /3-relaxation strength is 
difficult to understand. In particular, an explanation will 
be problematic if the IS dynamics in different regions of 
the PEL are similar as suggested by our results. 

In view of these findings, we propose that the MB 
transitions are the elementary steps of the a-process or, 
equivalently, the structural relaxation results from a se- 
ries of MB transitions. Further, we suggest that the /3- 
process is related to the exploration of a MB, i.e., to the 
E-sequences. In this picture, both relaxation phenomena 
are naturally multi-step processes. Further, the growing 
timescale separation of the a- and the /9-process upon 
cooling is a consequence of the increasingly important 
trapping of the system in a MB. When decreasing the 
temperature the system is forced into deeper and deeper 
regions of the MB so that smaller and smaller regions 
of the MB can be explored due to the /3-process and, 
consequently, its relaxation strength decreases. In other 
words, the number of IS visited during the /3-process is 
temperature dependent. Finally, the relative amplitude 
of the a- and the /3-process depends on the steepness of 
the MB and the number of IS combined in the MB. Thus, 



the experimental observation that the relative strength of 
both relaxations is different for various glass formers can 
be met. Of course, our speculations need to be care- 
fully checked in future investigations. As pointed out by 
Debenedetti and Stillinger^ , suitable data is not available 
at the present time. 

In conclusion, the following picture appears to emerge 
for IS and MB dynamics at sufficiently low temperatures. 
When the system jumps between IS organized into a MB, 
groups of spatially correlated particles achieve compara- 
tively large displacements by performing string-like mo- 
tion. In doing so, the strings arise, for the most part, 
not due to a coherent motion of all involved particles 
during a single IS transition, but from the interplay of 
displacements of various sub-units taking place at differ- 
ent times. During the exploration of the MB, the system 
enters energetically less favorable regions near the ridge 
of the MB from which it either escapes to a new MB 
or returns to the bottom of the old MB. While in the 
latter case, the formed strings dissolve due to the back- 
ward motions of the involved particles, they persist in 
the former. The particle rearrangements during different 
MB transitions are basically uncorrelated. In particular, 
different groups of particles are mobile and form strings. 
Hence, one may speculate that successful string-like mo- 
tion results in the escape from a MB. We ascribe the 
exploration of the MB to the Johari-Goldstein /3-process 
and transitions between distinct MB to elementary steps 
of the a-process. In this picture, both relaxation phe- 
nomena are multi-step processes and their growing time- 
scale separation on cooling results from the increasingly 
important temporary trapping of the system in the MB. 
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